Do small swimmers mix the ocean? 
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In this communication we address some hydrodynamic aspects of recently proposed drift mechanism of 
biogenic mixing in the ocean [K. Katija and J. O. Dabiri, Nature 460, 624 (2009)]. The relevance of the 
locomotion gait at various spatial scales with respect to the drift is discussed. A hydrodynamic scenario of the 
drift based on unsteady inertial propulsion, typical for most small marine organisms, is proposed. We estimate 
its effectiveness by taking into account interaction of a swimmer with the turbulent marine environment. Simple 
scaling arguments are derived to estimate the comparative role of drift-powered mixing with respect to the 
turbulence. The analysis indicates substantial biomixing effected by relatively small but numerous drifters, such 
as krill or jellyfish. 
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An intriguing question vigorously discussed in recent lit- 
erature is whether marine organisms contribute significantly 
to the ocean mixing. The arguments both in favor JHQ] and 
against |3] a substantial role of biomixing have been based 
on the overall energy balance, and assumed that a certain part 
of the overall energy intake of marine organisms contributes 
to mixing rather than being dissipated by other means. This 
approach may be inadequate in view of uncertainties in esti- 
mating alternative dissipation pathways without considering 
swimming dynamics explicitly. Very recent work l3l sug- 
gests that the net displacement of fluid particles by swim- 
mers (known as "drift mechanism" |5]), being enhanced by 
viscosity, is more relevant for mixing than biologically gen- 
erated turbulence. Experimental observation of drift behind 
swimming jellyfish and hydrodynamic demonstration estimat- 
ing drift capacity of passively towed objects at finite Reynolds 
numbers. Re, |4] give support to the notion of strong biomix- 
ing. It is the aim of this communication (i) to examine the hy- 
drodynamics of the drift induced by self-propelled objects; (ii) 
to derive mechanistic estimates based on scaling arguments, 
taking into account the turbulent environment and comparing 
swimmer's contribution to mixing with turbulent stirring on 
relevant scales. 

The major factors determining the answer to the question in 
the title are the drift volume D dragged by the swimmer and 
the turbulent scales characteristic of the marine environment. 
Both factors may vary within a wide range. 

A steadily towed body in a viscous fluid drags an infinite 
volume of fluid when the flow is dominated by viscosity |@], 
while in inviscid potential flow it is of the same order of mag- 
nitude as the volume of the object [5]. Recent demonstration 
of the drift due to a passively towed two-dimensional body 
of the dimension lb at Re = Uihjv ^ 5-100 in |4] in- 
dicates that the viscosity-enhanced drift volume may far ex- 
ceed the swimmer's own volume. This tendency can be read- 
ily explained. Advection of vorticity, shed by a body towed 
steadily with the speed V through a fluid quiescent at infin- 
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ity, is balanced at some distance downstream by transverse 
viscous diffusion. The asymptotic solution for the velocity 
deficiency in the wake v — U ~u{vs\ the coordinate frame co- 
moving with the swimmer in the negative x direction) can be 
derived assuming small departure from the free stream veloc- 
ity at positions sufficiently far downstream, t; <C C/, and us- 
ing the boundary layer approximation, djdx ^ d/dr, where 
r is the transverse radial distance. The asymptotic solution 
at a; — oo is i; ^ [J/^ttvx) exp (— [/r^/4i^a;) |7]. Here 
J5{x) is the point force yielding the kinematic momentum 
flux (force per unit mass) J = — Je, where e is the unit vec- 
tor along the x-axis. The flow disturbance outside the wake 
away from the object is UTOtational radial flow due to "poten- 
tial source", u ~ {J/A^^U)x/\x\^^, compensating for the mass 
deficiency in the wake fr]. The rate at which the fluid vol- 
ume is dragged across a stationary plane behind the swimmer, 
q = 2it vrdr = J/U, can be related to the net drag force 
exerted on the body, F = pj, where p is the fluid density 
Clearly, since J — const, the corresponding drift volume D is 
infinite, similar to the zero Re case |6]. Considering the par- 
tial drift volume Dp dragged by a body that had been towed 
a distance L over time L/U H, and expressing the drag as 
F — ^CdP U^A, where A is the wetted surface area and Cd 
is the drag coefficient, we find that Dp = qL/U = \ CdAL. 
Since Cd is a decreasing function of Re, such that Cd ^ 
Re~" with the exponent a = 1 at Re< 1, and varying from 
1/2 to 1/5, respectively, in the laminar and turbulent bound- 
ary layer regime at high Re, the concept of viscosity-enhanced 
drift behind passively towed body demonstrated numerically 
in [4] becomes clear, as lower Re results in the increased par- 
tial drift volume. The question is: to what extent this notion is 
relevant towards drift induced by self-propelled organisms? 

The flow field around self-propelled objects could be, how- 
ever, quite different from that of a passively towed body. A 
steadily self-propelled swimmer generates no net momentum 
flux, since the thrust is counter-balanced by the drag force 
JUlit]. The two forces of equal magnitude, acting in oppo- 
site directions and separated by some distance, constitute a 
force dipole, that can be approximated away from the body as 
Qij\7j6{x) (@]. Higher order terms, such as force quadrupole, 
MijkVj \7kS{x), etc., are decaying faster. The net momentum 
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flux induced by such force dipole is zero and therefore the cor- 
responding drift volume is expected to be far smafler than for 
a passively towed body. 

The drift induced by a high-Re steadily propelled swimmer 
can be estimated along the same lines as that by a passively 
towed object considered above. The far-field velocity defi- 
ciency in the zero-momentum wake behind a localized force 
dipole in a spatially uniform stream of velocity U along the x- 
axis is given by v = {Q/iniyx'^) (1 — 77)6"'' at a; cxd, where 
77 = Ur'^ / Avx 1 9] . The far flow field outside the wake is due to 
a potential doublet. It = (-(3/47r;7|a;|^){e-3(e-a;) x/\x\^}, 
and is the same as the far field of a steadily moving object in 
potential flow [7]. It can be readily shown that the net fluid 
flux q across a stationary plane in the wake is zero. Whereas 
for a monopole forcing exerted on a passively towed object the 
radial potential flow outside the wake is compensated by a net 
influx in the wake, for a dipole forcing the potential doublet- 
like flow outside the wake causes no net fluid drift, and, as 
result, there is no net drift within the wake as well. Thus, the 
net drift volume will be comparable to the swimmer's volume 
as in potential flow. 

Although small swimmers (<1 mm wide) are mainly non- 
vertical migrators |10], it is instructive to examine the drift 
associated with their motion. In the viscous regime, the mean 
far-field flow (axisymmetric with respect to x-axis) around a 
force-free swimmer propelled with the mean velocity U in the 
negative x-direction through unbounded fluid quiescent at in- 
finity, can be written in the comoving frame as 
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where Q and M are, respectively, the strengths of the force 
dipole and quadrupole, and r, 6 are spherical coordinates with 
measured from the positive x axis. A swimmer with neg- 
ative Q is called a "puller", as for micro-organisms having 
their thrust-generating apparatus in front of the body (which 
dominates the drag), e.g. biflagellate algae such as Chlamy- 
domonas. An organism with positive Q is a "pusher", i.e. the 
thrust is generated behind the body, as in bacteria or sperma- 
tozoa. Proceeding along the same Unes as in Ref. Da], it is 
convenient to express the radius of the streamline tube ^0 far 
up or downstream from the swimmer (which are equal due to 
fore-and-aft symmetry), through t/j = — [/^q/2. Expanding 
for large ^0 and solving for r yields: 
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The drift volume D is defined as the fluid volume between 
an undisturbed material plane which is initially far in front of 
the swimmer, and the final position of the same surface dis- 
turbed by the swimmer moving off to infinity The partial 
drift volume can be defined for a material surface of a finite 
radius ^0 ||6|1j Dp — lix X^ad^Q, where X is the displace- 
ment of a fluid particle in the direction of translation. Taking 
the limit of Dp at ^0 00 results in the drift volume D. 
Using the integral mass conservation arguments the par- 
tial drift volume can be found as Dp = A — V, where V 



is the swimmer's volume and A is the volume between the 
laterally displaced streamline tube and its unperturbed posi- 
tion, A = TT J^' (r^ sin0^ — ^g) da;, with Xi and Xo being the 
coordinates of the streamline tube inlet (upstream) and outlet 
(downstream). Using r(^o) in Eq. (|2]l, we change the integra- 
tion variable to 6 via x = ^0 cot and evaluate the integral as 
A = AMt:/Uu + 0(^(7^) at 00. Thus, the drift volume 
for a force-free viscous swimmer is finite. 
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Note that the drift volume is independent of Q whereas both 
"pushers" and "pullers" would displace the same amount of 
fluid. The value of M and the corresponding drift volume 
depends on particular swimming gait. For instance, it can 
be evaluated for a model "spherical squirmer" ifTlll of ra- 
dius a propelled by purely tangential axisymmetric and time- 
independent surface motion, that in the co-moving reference 
frame can be written as M|,=a = Us(l + l3cosff)sm9eg, 
where Us is the magnitude of the surface velocity. The propul- 
sion speed and dipole and quadrupole forcing intensity in ([T]) 
can be readily found upon application of the latter boundary 
condition at r = a, yielding U — 2us/3, M = a^Ugv/i, 
Q = —a^(5usv/2. Therefore, the drift volume for a "spher- 
ical squirmer" is finite and equal to half of its volume, D = 
AMtt/Uv — V = V/2. The same result was found for a rigid 
sphere moving through an inviscid fluid |5]. 

Although the above analysis shows no considerable drift 
and the anticipated drift volume is comparable to that found 
for inviscid potential flow, yet a substantial drift observed in 
experiments with dye dispersion behind a jellyfish f7| can 
be attributed to unsteady locomotion, typical for small in- 
ertial marine swimmers lll2|] . For an organism that swims 
unsteadily, such as jellyfish propelled by jetting, or small 
crustaceans propelled by rowing their legs, the instantaneous 
thrust must be balanced by the force to overcome its inertia 
as well as two hydrodynamic forces; drag and the acceler- 
ation reaction (or the added mass force) lfl2l - [l4ll . During a 
fast propulsion stage (power stroke), the swimmer transfers 
momentum to the surrounding fluid generating accelerated 
propulsion, as the thrust exceeds the drag force and the reac- 
tive force by the amount equal to the swimmer's inertia. Dur- 
ing a slow retraction stage, the swimmer decelerates absorb- 
ing back some momentum, but the overall net momentum flux 
over a stroke cycle transferred to the fluid is non-zero due to 
a nonlinear nature of the inertial propulsion. This may yield a 
considerable drift due to a momentum/mass deficit in the wake 
in the same way as for a passively towed object discussed 
above. These arguments may shed light on why numerical 
demonstration of viscosity-enhanced drift behind a passively 
towed object, provided in |4] without a comment or justifica- 
tion, is relevant, to some extent, to marine self-locomotion. 

In turbulent environment, the drift volume is not infinite, 
as it is dispersed by the background turbulence. The turbulent 
marine environment is characterized by the Kolmogorov scale 
4 = [v^/eY/'^, where V is the kinematic viscosity and e is 
the dissipation rate per unit mass of the fluid, and the Ozmidov 
scale, to = (e/iV'^)^/^, corresponding to the largest eddy that 
can overturn for prescribed stable density stratification in the 



3 



ocean dp/dz. Here N = -y/— (g/po)(dp/dz) is the buoyancy 
frequency, with po being the mean fluid density [15]. 

When the swimmer is propelled though a weakly turbulent 
fluid, the momentum wake fades away at a distance £ down- 
stream where the velocity deficiency v becomes comparable 
to the magnitude of turbulent velocity fluctuation u' . More- 
over, the vorticity in the wake is diffused transversely by tur- 
bulent eddies with the the scale-dependent turbulent disper- 
sion coefficient, JC. Although the detailed analysis of the wake 
in unsteady locomotion in a weakly turbulent environment is 
beyond the scope of the present study, a rough estimate of the 
mean momentum flux yields J = {m/ p){AU/ At), where 
AU is the velocity gain by a swimmer of mass m during a 
pulse phase of duration At. Taking constant eddy diffusivity 
corresponding to the largest eddy, and assuming that the time- 
averaged structure of the momentum wake is the same as for 
a passively towed body, the comparison between the veloc- 
ity deficit in the momentum wake responsible for the drift, 
V ~ J/K.£, and u' yields the following scaling estimate for 
the apparent wake length, I, 

e - {mAU)/{ATpJCu') . (4) 

For instance, the jellyfish Aurelia aurita with the bell diam- 
eter of ^ 8.5 cm and mass m ~ 20 gr, gains AU ^ 3 
cm s~^ during At ~ 0.5 s lITill . For a typical turbulent 
patch in a seasonal thermocline with e = 10^^ W kg^^ and 
buoyancy frequency N ^ 0.01 s~^ lfl6ll . the turbulent ve- 
locity fluctuations at the Ozmidov scale, u'^ ~ {eio)^^^ are 
about 1 mm s^^, and the corresponding vertical eddy diffu- 
sivity (also known as density or diapycnal diffusivity, fistl) 
K. - u'Jo = e/N^ ~ 10-* m2 s^^. Thus, the estimate dUl 
yields a surprising £ ^ 10 m! However, the strength of the 
momentum flux diminishes rapidly with the organism dimen- 
sions as expected. For a smaller jellyfish with the bell diame- 
ter of 1.8 cm and mass 0.5 gr, which gains AU ^ 1.2 cm s^^ 
in At ^ 0.3 s I14<1 . the estimate (|4|i yields ^ ^ 20 cm only. 

Note that the proposed mechanism leading to (01 does not 
offer any viscosity-enhancement of the drift. On the opposite, 
smaller swimmers are expected to be less efficient drifters. 
For small viscous swimmers (propelled either steadily or un- 
steadily) the proposed mechanism is not operative. However, 
tiny swimmers can be better drifters when they swim against 
an external force, such as gravity. In this case they generate 
a slowly decaying Stokeslet in the far field, analogously to 
a passively towed object |6]. This situation arises e.g. dur- 
ing upward active swimming or downward passive sinking of 
bottom-heavy swimmers. The net momentum transferred to 
the fluid will be proportional in this case to the swimmer's ex- 
cess weight and will be aligned with the direction of gravity 
regardless of the direction of their motion. 

The volume of a typical size £ > £k is dragged along by a 
swimmer, which can be considered moving independently of 
the weakly turbulent local flow [161 , for a time r which takes 
for turbulent pulsations on the same scale £ to disperse it. The 
estimate for this time is t ~ £/u', where u' is the charac- 
teristic turbulent fluctuation velocity on the scale £. It can be 
estimated, assuming it obeys the isotopic Kolmogorov scaling 
in the inertial sub-range, as u' ^ (e^)^/-^. The resulting effec- 



tive dispersion coefficient is 2? U'^t ^ (j) U'^£/u', where 
U is the swimmer velocity and 4> is the volume fraction occu- 
pied by the swimmers. This can be compared with the vertical 
eddy diffusivity K, ^ s/N'^, leading to an order-of-magnitude 
estimate 

V/K. - (/) U^N^£^/^e-'^^^ . (5) 

For small viscous swimmers (Re= U£s/v < 1) of a size 
^ ^k, with a swimming speed lower than the rms turbu- 
lence velocity, such as phytoplankton lfT6l[ . the typical size of 
the dragged volume accessible to dispersion by turbulent pul- 
sations is comparable to £k. The typical dispersion time is 
defined in this case on the finest turbulence scale and is esti- 
mated as T - [v/ef/'^ leading to P - U'^{v/ef/'^. This 
should be compared with the turbulent dispersion coefficient 
JC to give 

Note that the latter estimate does not depend on the swimmer's 
size. 

For large inertial swimmers corresponding to Re ^ 1, 
where the potential flow assumption is applicable, £ ^ £s 
[5] and the dispersion time is restricted to r ^ ^s/C^, as a 
tracer particle interacts with a swimmer for a typical time that 
takes a swimmer to pass by. This modifies the dispersion es- 
timate XoT) ^ (f) U£s. This scaling-based estimate coincides, 
up to an 0(1) shape-dependent multiplicative constant, with 
the results [flST of rigorous analysis of stirring by 2D and 3D 
bodies propelled in the potential flow, which have shown that 
displacement of passive tracer particles is completely domi- 
nated by "head-on" collisions with the bodies. Comparing to 
vertical dispersion by turbulent buoyancy flux yields 

V/K.^ (j)U £sN^e^^ . (7) 

However, the assumption of potential flow does not take into 
account a wake shed by the self-propelled body and the pre- 
diction Q may considerably underestimate the contribution 
of drift-induced mixing. 

The above estimate of turbulent dispersion /C ^ e/N"^ sug- 
gested by scaling is valid in the transitional regime corre- 
sponding to strong density stratification, where the parameter 
/ = e/vN"^ {i.e. the squared ratio of the buoyancy time scale 
N^^ to the time (i//e)^/^ required for turbulent events to de- 
velop) is in the range 7 < / < 100 [17]. Beyond this range 
JC ^ e/TV^ greatly overestimates the actual turbulent disper- 
sion. In the energetic regime corresponding to weak strati- 
fication, where / > 100, turbulence is only weakly affected 
by buoyancy as it develops relatively fast, and the dispersion 
coefficient is /C - v{e/vN'^ f''^ Oj. The use of the latter 
scaling estimate for JC modifies the ratio V/JC in Eqs. (|5])-(|7l) 
resulting in even stronger contribution of the drift to mixing 
in the energetic regime. 

Let us estimate the relative effect of mixing caused by verti- 
cally migrating Antarctic krill Euphausia pacifica. The buoy- 
ant frequency N for a strongly stratified ocean surface layer is 
typically about 10^^ s^^ and the background dissipation rate 
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measured during daylight when krill remained stationary at 
depths - 100 m is £ « IQ-s W kg-^ 1 1]. This gives / ~ 10, 
and therefore ^ is applicable. As the detailed measurements 
of propulsion kinematics available for jellyfish (e.g. fl?]) 
have not been reported for krill, we cannot estimate £ using 
(|4|l. We, therefore, take £ w £s as a conservative estimate of 
the characteristic dragged volume dimension. The krill length 
£s ~ 1.5-2 cm, the average swimming speed is [/ w 2.5-3.5 
cm s^^, and volume concentration for a dense layer of krill 
is = 0.07-0.14%, con-esponding to 5, 000-10, 000 individ- 
uals m~^ fit], whereas the volume occupied by a single krill 
is estimated as Trr^^^, with the krill radius is r « 0.15 cm 
ll20ll . This yields V/K. ^ 5-15 in Eq. Q, which indicates 
a substantial biogenic contribution of the drift to the overall 
mixing. 

In the seasonal thermocUne stratification the buoyant fre- 
quency N may exceed 0.05 s^^ at the depth between 40 and 
240 m so that turbulence is effectively suppressed, and 
biomixing is expected to be a major stirring mechanism. In 
this case, / ~ 0.4 corresponds to the mo/ecM/ar regime, when 
IC should be replaced in the comparison estimates by f ITtIi . 
In the latter case, V ^ (p U'^t, where r ~ ^"^l^- There- 
fore, V/v ~ U'^£^ jv^. For Antarctic krill [1], we have 
T)lv ~ 50-250. On the other hand, for the weakly strati- 
fied ocean, N ^ Vir^ s^^ and less, and assuming e ~ 10~^ 
W kg"^ we obtain / ~ 1000. This suggests that the ratio 
D jK should be estimated using K ~ v{e/vN'^Y/'^ corre- 
sponding to the energetic regime and yielding V/K. ^ 2-4. 

Most of the marine biomass is concentrated in small organ- 
isms, such as copepods (£s ~ 1-3 mm), with some of them ex- 
hibiting vertical migration, feeding near the surface at night, 
then sinking into deeper water during the day to avoid visual 
predators 1 10]. The question of whether such diel vertical mi- 
grations may contribute to significant dispersion enhancement 
is still under discussion (e.g. ]^). According to the estimates 
in Eq. (|6]l, small swimmers of the size £s < £k may gener- 



ate significant mixing only at relatively high densities, due to 
their small propulsion speed. Using typical copepod swim- 
ming speeds u ~ 2-5 mm s^^ at a maximum volume con- 
centration 1.5-2% (based on the universal scaling for maxi- 
mum packing in a swarm ] 19]) and the same values of e and 
N as in Ref. |1], Eq. (|6j yields V/K. ~ 1 in the transitional 
regime. In the energetic regime with N = 10~^ s~^, Eq. (|6|l 
gives V/K ^ 0.1. The mean population density of copepods 
in the ocean is estimated as n < 10"* m~^ 112 ill , correspond- 
ing to (j) < 0.001%. For such low mean population density, 
the above estimates suggest that drift-induced contribution to 
mixing due to copepods vertical migration is negligible. 

The drift may be, however^facilitated by potential collec- 
tive effects, as suggested in The drift volumes of indi- 
vidual swimmers can be shared by many individuals within a 
swarm, so that coordinated migration may be responsible for 
the formation of a large collective trail on the swarm scale. 
Assuming regular spatial orientation of the drifters, one can 
estimate the mean nearest-neighbor distance in a swarm as 
~ 1.2 n~^/^ ll22ll . Using the universal scaling for maximum 
population density fl3l, rim ^ 0.022 £~^'^ one arrives at the 
mean nearest-neighbor distance in a dense swarm ^ 4.25 £\''^. 
This estimate suggests that organisms with the body length '-^l 
mm in a dense swarm are separated by about one body length 
on the average, while for organisms with £s ^ ^ cm, the mean 
nearest-neighbor distance is about two body lengths. These 
estimates are in accord with population densities reported in 
the literature and calculated from photographic recording of 
aggregations, acoustic measurements and catches in nets and 
trawls. For the organisms in the size range £s ~ 1-3 cm, 
such as Euphausia pacifica and others, the observed nearest- 
neighbor distance in dense swarms is in the range of 2— 4 body 
lengths (see Tab. 2 in ]23]). 

The intriguing problem of collective drift by hydrodynam- 
ically interacting self-propelled objects will be subject of fur- 
ther investigation. 
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